OBSERVATION OF SHELL STRUCTURE IN SODIUM NANOWIRES 
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The quantum states of a system of particles in a finite spatial domain in general con- 
sist of a set of discrete energy eigenvalues; these are usually grouped into bunches of 
degenerate or close-lying levels pj, called shells. In fermionic systems, this gives rise to 
a local minimum in the total energy when all the states of a given shell are occupied. 
In particular, the closed-shell electronic configuration of the noble gases produces their 
exceptional stability. Shell effects have previously been observed for protons and neu- 
trons in nuclei and for clusters of metal atoms [w-W. Here we report the observation of 
shell effects in an open system — a sodium metal nanowire connecting two bulk sodium 
metal electrodes, w^hich are progressively pulled apart. We measure oscillations in the 
statistical distribution of conductance values, for contact cross-sections containing up 
to a hundred atoms or more. The period follows the law expected for the electronic 
shell-closure effects, similar to the abundance peaks at 'magic numbers' of atoms in 
metal clusters rtaJ4l. 
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Metallic constrictions, in the form of nanowires con- 
necting two bulk metal electrodes, have recently been 
studied down to sizes of a single atom in cross section by 
means of scanning tunnelling microscopy (STM) and me- 
chanically controllable break-junctions (MCB) g,^. By 
indenting one electrode into another and then separat- 
ing them, a stepwise decrease in electrical conductance 
is observed, down to the breakpoint when arriving at 
the last atom. Each scan of the dependence of conduc- 
tance G on the elongation d is individual in detail, as the 
atomic configuration of each contact may be widely dif- 
ferent. However, statistically, many scans together pro- 
duce a histogram of the probability for observing a given 
conductance value, which is quite reproducible for a given 
metal and for fixed experimental parameters. 

A simple description of the electronic properties of 
metallic nanowires is expected to work best for mono- 
valent free-electron-like metals. The alkali metals are 
most suitable, as the bulk electronic states are very 
well described by free particles in an isotropic homo- 
geneous positive background. In a previous experiment 
on sodium point contacts W a histogram was obtained, 
which showed pronounced peaks near 1, 3, 5 and 6 times 
the quantum unit of conductance, Gq = 2e^ /h (where e 
is the charge on an electron and h is Planck's constant). 
This is exactly the series of quantum numbers expected 
for the conductance through an ideal cylindrical conduc- 
tor. We have now extended these experiments to higher 
temperatures and a much wider range of conductance 
values. 

The MCB technique |^ was adapted for the study of 
nanowires of the highly reactive alkali metals following 
Ref. (see Fig. 1). Scans were taken continuously by 
ramping the displacement d of the electrodes with respect 



to each other, using the piezo-electric driver. Each indi- 
vidual curve of conductance versus displacement, G{d), 
was recorded in ~0.1 seconds from the highest conduc- 
tance into the tunnelling regime. Histograms of conduc- 
tance values were accumulated automatically involving 
lO'^-lO^ individual scans, and for a number of sample 
temperatures between 4.2 K and 100 K. Reproducibility 
and reversibility were verified by periodically returning 
to the same measuring conditions. Here we present re- 
sults for sodium, but similar results were obtained for 
lithium and potassium. 

Figure 2 shows the temperature dependence of the his- 
tograms for sodium in the range from to 20 Gq. In order 
to compare different graphs, the amplitude has been nor- 
malised by the area under each graph. The histograms 
at low temperatures and low conductances are similar to 
those given in ref. m. We clearly recognise the familiar 
sharp peaks below 6 Go, while at higher conductances a 
number of rather wide maxima are found. With increas- 
ing temperature we observe a gradual decrease in am- 
plitude of the lower conductance peaks, while the high- 
conductance peaks dramatically sharpen and increase in 
amplitude. 

Our central results are shown in Fig. 3. The histogram 
for sodium at temperature T = 80 K up to G/Gq = 
120 is shown in Fig. 3a inset. The smooth background 
shown by the dotted curve was subtracted, in order to 
present only the oscillating part; the latter is plotted on 
a semi-logarithmic scale in the main panel of Fig. 3a. 
Up to 17 oscillations are observed, at positions which are 
reproduced well for ^ 10 different samples. The influence 
of the procedure of background subtraction on the peak 
positions is much smaller than the width of the peaks. 

The radius. R, of the narrowest cross section of the 



nanowire can be obtained from the semi-classical expres- 
sion for a ballistic wire [pi , 
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where kp is the Fermi wavevector. Using this expression 
we calculate the radius of the nanowire at the peak po- 
sitions (in units of kp^)] the results are shown in Fig. 3b 
(open squares) for consecutively numbered peak index 
i. Thus we observe that the peak positions are peri- 
odic as a function the radius of the wire. This period- 
icity suggests a comparison of our experimental results 
with the magic numbers observed for sodium clusters. 
In the cluster experiments, a remarkable structure was 
observed in the distribution of cluster sizes produced in 
a supersonic expansion of sodium metal vapour p,S-n3| . 
In the spectrum of the abundance versus the number of 
atoms, N ^ in the cluster, distinct maxima are observed 
for A^ = 2, 8, 20, 40, 58, 92, .... These 'magic numbers' cor- 
respond to clusters having a number of valence electrons 
(equal to the number of atoms) which just complete a 
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FIG. 2. Temperature evolution of sodium histograms in 
the range from to 20 Gq. The voltage bias was 10 mV 
and each histogram is constructed from 1000-2000 individual 
scans. The amplitude has been normalised by the total area 
under each histogram (a.u., arbitrary units). 
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FIG. 1. Schematic view of the MCB technique for alkali 
metals |7|. While immersed in paraffin oil, the sample is 
pressed onto four 1 mm diameter brass bolts, which are glued 
onto the isolated substrate. Current and voltage leads are 
fixed with nuts onto the bolts. A notch is cut into the sam- 
ple at the centre. This assembly is mounted inside a vacuum 
can, which is immersed in liquid helium. By bending the 
substrate at 4.2 K in vacuum, the sample is broken at the 
notch. Two fresh surfaces are exposed, and their distance is 
adjusted by changing the force on the substrate, employing a 
piezo element for fine control. The temperature of the sample 
is controlled by a heater and thermometer, which are fixed at 
the bottom of the substrate. The cryo-pumping action of the 
low-temperature environment ensures that the surfaces are 
not polluted by adsorbates. The purity of the sodium metal 
was at least 99.9%. The conductance was recorded using a 
standard d.c. voltage bias four-terminal technique, measuring 
the current with 16-bit resolution. The drift and calibration 
of the current-to-voltage converter was verified against stan- 
dard resistors corresponding to 1, 10 and 100 Go, leading to 
an overall accuracy in the conductance better than 1% for 
G > lOGo. 



shell, where the wave functions arc considered as freely 
propagating waves inside a spherical potential well. The 
magic numbers are periodically distributed as a function 
of the radius of the clusters, where the radius is obtained 
from the number of atoms, N, in the cluster using [g|j] 
kpR = 1.92A^^/'^. In Fig. 3b we compare the radii for the 
cluster magic numbers (filled circles) with the radii cor- 
responding to the peaks in the conductance histogram 
(open squares). We establish a direct correspondence 
between preferential values for the conductance in his- 
tograms for sodium nanowires, and the cluster magic 
numbers. 

The linear relation between the radii at cluster magic 
numbers and the shell number arises owing to fluctua- 
tions in the density of states as a function of kpR, which 
in turn give rise to fluctuations in the free energy of the 
system, J7(i?). The periods of oscillation can be expressed 
in terms of semiclassical closed orbits inscribed inside the 
boundaries of the system ||l|,Q| , where the path length L 
should be an integer multiple of the Fermi wavelength 
Xp = 2'K/kp. Fluctuations in the density of states for a 
free electron nanowire and their influence on the free en- 
ergy has been recently considered in a number of theoret- 
ical papers ||lj-|l^. Using the semiclassical expansion of 
the free energy for a cylindrical wire [|l8|,|9| including only 
the three lowest order terms, we calculated the positions 
of the energy minima, which are plotted as a function of 
shell number in Fig. 3b (open triangles). The agreement 
with the experimental points is fairly good, and may be 
improved taking into account a more realistic shape of 
the potential well and higher-order contributions. 

Our interpretation uses the approximately linear semi- 
classical relation G{R) between conductance and wire 



cross-section, equation (1). Corrections may arise from 
two mechanisms. First, back scattering on defects (or 
phonons) may shift the peak positions. The fact that 
strong scattering would tend to smear the peak struc- 
ture, and the close agreement between the experimental 
and theoretical periodicity both suggest that this shift is 
small. Second, the mechanism of level-bunching leading 
to fluctuations Q{R) should also give rise to fluctuation 
corrections to G{R), which by itself would lead to peaks 
in the histograms even for a perfectly smooth distribu- 
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FIG. 3. Conductance histograms for sodium, showing ev- 
idence for shell structure, (a), Histogram of the number of 
times each conductance is observed versus the conductance in 
units of Go ~ 2e^ /h. These data are for sodium at T = 80 K 
and bias voltage V — 100 mV, constructed from over 10,000 
individual scans. Inset, the raw data and the smooth back- 
ground (dotted curve); the background is subtracted to give 
the curve in the main graph. The logarithmic scale for ordi- 
nate axis helps to display the smaller-amplitude features at 
high conductance values, (b), Radius of the nanowire at the 
positions of the maxima in (a) versus peak numbers (open 
squares), where R is given in units fcp^). The radii at the 
peak positions are compared with the radii corresponding to 
the magic numbers for sodium metal clusters (filled circles) 
Jill and with those expected from a semiclassical description 
for the fluctuations in the free energy for the nanowire (open 
triangles) OM. 



tion of wire diameters. Peaks due to this mechanism 
would be found at those points where the increase of the 
conductance with wire diameter is slow, and indeed we 
believe the peaks at low conductance, near 1, 3 and 6 Go 
are due to this mechanism. However, we argue that the 
main structure in Fig. 3 is due to the shell structure in 
n{R). 

The main argument in favour of this interpretation 
comes from the temperature dependence shown in Fig. 2. 
As in the cluster experiments pi, the electronic shell 
structure becomes observable at higher temperatures, 
where the increased mobility of the atoms allows the 
system to explore a wider range of neighbouring atomic 
configurations in order to find the local energy minima. 
At low temperatures, the nanowires are frozen into the 
shape, which evolves from mechanical deformations in 
the breaking and indenting processes. The decrease in 
amplitude of the sharp quantization peaks at low conduc- 
tance can probably be explained by thermally induced 
breaking of the contact when it consists of only a few 
atoms. If the fluctuations in G{R) were solely responsi- 
ble for the observed peak structure, it could be argued 
that the higher temperature would favour formation of a 
smoother and longer wire, leading to less backscattering 
of electrons and less tunnelling corrections, respectively. 
Backscattering is responsible for the shift to lower val- 
ues of the conductance peaks near 1, 3, 5 and 6 Go [ pO[ , 
and we find that the position of these peaks does not 
change with temperature. An indication of the length 
of the wire may be obtained from the global variation of 
the conductance with elongation pl| and the time evo- 
lution for fixed elongation p2| . In our experiment, both 
variations show that the effective wire length decreases 
for increasing temperature. With these arguments an in- 
terpretation of the temperature dependence of Fig. 2 in 
terms of G{R) fluctuations alone is ruled out. 

The correspondence between the shell structure in clus- 
ters and in nanowires can be explained by comparing the 
energy levels for a three-dimensional spherical potential 
well and a two-dimensional cylindrical geometry. The 
levels are obtained from the zeros of the spherical and or- 
dinary Bessel functions, respectively. Apart from a small 
constant shift, the distribution of levels for the two sys- 
tems has a very similar structure, with gaps and bunches 
of levels occurring at the same positions. This explains 
the striking similarity of the magic radii for the clusters 
and nanowires observed in Fig. 3b. 
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